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We consider a charged spinless quantum particle confmed to a graph 
consisting of a loop to which a halfline lead is attached; this system 
is placed into a homogeneous magnètic field perpendicular to the loop 
plane. We derive the reflection amplitude and show that there is an 
infinite ladder of resonances; analyzing the resonance pole trajectòries 
we show that half of them turn into true embedded eigenvalues provided 
the flux through the loop is an integer or halfinteger múltiple of the flux 
unit hc/e. We also describe a general method to solve the scattering 
problem on graphs of which the present model is a simple particular case. 
Finally, we discuss ways in which a state localized initially at the loop 
decays. 



I begun my career at times when the world was much less connected, and of most 
parts we knew only from journals arriving not quite regularly. My first encounter 
with the mathematical theory of unstable quantum systems and related scattering 



problems occurred in this way, particularly through [24. 25] and related papers by 
Larry Horwitz; the subject remained for me as well as for him an old love to which 
we return regularly from time to time. Only many years later I had an opportunity 
to met him in person and to appreciate also his charisma. A distinguished birthday 
is a good opportunity to come with another decay scattering system; I should say 
that I prefer presents which are amusing rather than expensive. 



1 Introduction 

Quantum mechanics for a nonrelativistic particle whose configuration space is a 
graph has been studied already a long time ago in connection with models of orgànic 
molècules |33|] . A recent new interest to these problems 0, ü, |, H, 0, |Ï2|, [R| [L7], |20 



23| has been stimulated, in particular, by the progress of experimental sòlid state 
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Figure 1: A lasso graph 



physics which allows us to produce semiconductor "quantum wire" structures and 
other "mesoscopic devices" ; quantum mechanical graphs represent a natural model 
for many of them. 

Graph systems provide a convenient mean to study various quantum effects both 
from the theoretical and experimental points of view, because the freedom in setting 
the geometry of the configuration space allows us to create different dynamics; at 
the same time, models of this type are often explicitly solvable. This concerns, 
in particular, resonance scattering effects associated with the existence of quasi- 
stationary states in graph loops and appendices — see, e.g., JL8|. These effects fit 
well, of course, into the general theory of decay scattering systems as exposed in 
TÜ| , Chaps. 1,3], HH, or [|32], Sec. XII. 6], but they also make it possible to illustrate 



it and to draw fully specific conclusions. 

Our aim here is to investigate one more solvable model of this type. It consists 
of a halfline attached to a loop placed into a magnètic field; the parameters are the 
magnètic flux through the loop and "coupling strengths" between the graph links 
at the junction. Our analysis differs from an earlier treatment of similar systems 
II [ÏÏJ in several aspects. First of all, we consider a different and more general 
coupling between the loop and the halfline, and we put emphasis on the analytical 
solution of the problem. Furthermore, we shall be concerned with the decay and 
scattering properties of the system rather than with persistent currents induced by 
the magnètic field. 

Let us review briefly the contents of the paper. The Hamiltonian of the model 
we are going to study is introduced in the following section. Next, we derive in 
Section 3 its spectral and scattering properties. Then we make a digression and 
describe a general method to treat scattering problem on an arbitrary graph by 
"discretizing" it, í.e., transforming the corresponding Schrodinger equation into a 
set of linear equations involving just the wavefunction vàlues at the graph nodes. 
Returning to our model, we analyze in Section 5 its resonance structure by deriving 
an explicit expression for the resolvent and finding the resonance-pole trajectòries. 
Finally, in the concluding section we treat our model as a decay system and show 
how a state localized initially at the loop decays (or does not decay) in the course 
of time. 



2 Description of the model 

Consider a quantum particle confined to the lasso-shaped graph T sketched on 
Figure 1, i.e., a circular loop of radius R to which a halfline lead is attached. We 
suppose that the particle is nonrelativistic, spinless, and charged. To be specific 
we assume that its charge is q = — 1 ; we adopt the usual rational system of units, 
e = c = 2m = h = 1 . The system is placed into a homogeneous magnètic field 
of intensity B ; the vector potential A can be chosen tangent to the loop with the 
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modulus 

1 $ 

A =2 BR= L' (2 - 1} 
where $ is the magnètic flux through the loop and L is the loop perimeter. Under 
the convention we have adopted the natural flux unit || is hc/ e = 2tc , so the rhs 
of fl2.1|) can be also written as (p/R where <fi is the flux value in this scale. 

The state Hilbert space of the model is U = L 2 (T) := L 2 (0, L) © L 2 (1R + ) ; the 
wave functions will be written as columns, ip — (j \ . To construct the Hamiltonian 
we begin with the operator describing the free motion on the loop and the lead under 
the condition that the graph vèrtex is "fully disconnected" , so = Hi oop (B) © 
#haiflinc , where 

H loop (B) = (-td x + A) 2 , # halfline = -d 2 x (2.2) 

with the Dirichlet condition, u(0) = u(L) = /(0) = ; if there is no danger of 
misunderstanding we abuse the notation and employ the same symbol for the arc- 
length variable on both parts of the graph. The operator iíi oop has a simple discrete 
spectrum; the eigenfunctions 

Xn ^ = VnÏÏ V2£7 ' n = 1 i 2 T-- (2-3) 
correspond to the eigenvalues (^) , which are embedded into the continuous spec- 



trum of -f/haifline covering the interval [0, oo) . Notice that the effect of the magnètic 
field on the disconnected loop amounts to a unitary equivalence, 

#ioop(5) = U„ A H loop (0)U A , (2.4) 

where (Uau)(x) := e lAx u{x) . 

To couple the graph parts one has to follow the Standard strategy [0 which 
means to replace Dirichlet by a "connected" boundary condition at the vèrtex. In 
general, there is a nine-parameter family of such conditions. This is too many; we 
will be concerned with its three-parameter subfamily [17], |I8| , in particular, with 
a one-parameter set of boundary conditions known as 5-coupling [0. Hence the 
Hamiltonian of our model acts as the free operator specified by 



H^B)§ = (-^ A ;; +A ' U )-, (2.5) 

the wave function is continuous on the loop, 

u(0) = u{L) , (2.6) 

and satisfies the requirements 

/(0) = um(0)+/i/'(0), 



(2.7) 



u'(0)-u'(L) = af(0)-uf'(0) 
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for an a, [i G 1R and ui e (C; the vàlues of w, / and their derivatives at the vèrtex are 
understood as the appropriate one-sided limits. However, we shall restrict ourselves 
to the case of time-reversal invariant couplings which means to assume that to is 
also real; it has the meaning of a coupling constant between the loop (with a point 
interaction) and the halfline. In physical terms the conditions ( |2.6|) and ( |2.7| ) express 
the conservation of probability flow at the junction. 

The í-coupling corresponds to the choice [i = and cu = 1 in which case the 
wavefunction is fully continuous, 

u(0) = u{L) = /(0), (2.8) 

and 

«'(0)-u'(L)+/'(0) = a/(0); (2.9) 

for the sake of simplicity we shall write H a = H a) o,i • The parameter a is a coupling 
constant between the disconnected loop and the halfline; the fully decoupled case 
corresponds to a = oo as the notation suggests. 

Remarks. (a) The choice of the coupling at the vèrtex corresponds to a conceivable 
quantum-wire experiment. There is an approximation result [OI which means that 
the í-coupling constant a can be regarded as a mean value of a sharply localized 
potential. This corresponds, e.g., to a screened elèctrode placed at the vicinity of 
the junction; in a similar way one can model some of the more general boundary 
conditions ( |2.6| ) and ( |2.7| ) relating the parameters to physical quantities which an 
experimentalist can tune. 

(b) In general, the vector potential enters the boundary conditions — see || 
and the remarks in Section 4.3 below. In the present case, however, the outward 
tangent components of A at the junction have opposite signs, so their contributions 
cancel. This may not be true if the loop is noncircular and has còrners or cusps, 
but one can always achieve a cancellation by a suitable gauge choice. If the loop is 
viewed from outside as in the scattering process, the only quantity which matters is 
the magnètic flux <£> threading it. 

(c) The S-matrix for a coupling of three semiinfmite wires equivalent to Q2.6| ) 
and Q2.7| ) was derived in |fT?| . This comparison shows, in particular, that choosing 

\i = and putting e := (2+^2) > we obtain the coupling used in ||. On the 



other hand, the authors of [30] worked with the ideal 5-coupling, a = . 



3 Scattering and bound states 

Consider now the scattering problem on T , i.e., the reflection of a particle traveling 
along the halfline from the magnetic-loop end. We limit ourselves to the stationary 
formulation looking for generalized eigenvectors, in other words, solutions of the 
equation H a (B)ip = k 2 ip which satisfy the defmition domain requirements with 



exception of global square integrability. In view of ( |2.5| ), the most general Ansatz 
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for such a solution is 



uix \ 



p e~ iAx sin(kx + 7) 



/(*) 



e -ikx +re ikx 



(3.1) 



with ( fc-dependent) parameters r, (3 , and 7 ; the latter is generally complex. 

To find them we employ the boundary conditions. The identity ( |2.6|) in combi- 
nation with (2A) leads to the relation 



tan 7 



sin kL 



(3.2) 



The conditions 
solved by 



e 1 * — cos kL 

yield then a system of two linear equations for r, (3 which is 



(1 + ikfi) 


a — 


n ' 

sin 7 


+ iív 2 k 


(1 — ikfi) 


a — 


n ' 

sin 7 


— iuj 2 k 



with 



n 

sin 7 



k cos 7 — iA sin 7 — e * [k cos(kL + 7) — %A sin(kL + 7)] . 



Using again ( ^6|) and (|3.2|) , we arrive after a simple àlgebra at the expression 
r(k) = — 



(1 + ikjj) 


a — 


2k 
sin kL 


(cos $ - 


cos kL) 


+ ícj 2 A; 


(1 — ikfi) 


a — 


2k 
sin kL 


(cos $ - 


cos fcL) 


— iuj 2 k 



(3.3) 



for the reflection amplitude, in particular, 

{oí + ik) sin fcL — 2fc(cos $ — cos fcL) 



r(k) = — 



a — ik) sin kL — 2/c(cos $ — cos fcL) 



(3.4) 



in the 5-coupling case. This (lxl) S-matrix can be also written by means of the 
phase shift. For instance, denoting 



A(k) = A(a, <£>;&) := (a - ik) sin kL - 2fc(cos $ - cos kL) 

we can write the rhs of ( |3.4j ) as e 2í<5< - fc - ) with 

7r k sin kL 

5{k) = - +arctan^ J . 



(3.5) 



(3.6) 



As usual the growth of the phase shift is related to the number of scattering reso- 
nances within a given energy interval. It is clear from (|3.6| ) that 5(k) passes odd 
múltiples of 7r/2 whenever the denominator ( |3.5|) passes zero, of course, when there 
is not a simultaneous zero in the numerator. The last named situation occurs if 
and only if the flux $ through the loop is a múltiple of n . Hence "one half ' of 
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resonances is missing in that case; similar conclusions can be made in the general 
case of boundary conditions ( |3.3| ) when 



7T 



u 2 k 



S(k) = h arctan < ak + ^ — ; > . (3.7) 

2 l <* ~ ^k( cos ® ~ coskL ) J 

This is related to the existence of embedded eigenvalues at integer/halfinteger vàlues 
of <j) which will be clearly seen from the resonance pole trajectòries discussed below. 
The bound states can also be found directly: 

(a) It is clear that positive-energy bound states may be supported only at the loop. 

If we restrict our attention to the nontrivial case u 7^ , this is possible in 
view of (|2.7|) when u(0) = u'(0) — u'(L) = . Hence such bound states exist 
only at integer/halfinteger vàlues of the magnètic flux (in the natural units) 
and the corresponding eigenfunctions are given by (|2.3| ) with an even n for 
integer and odd n for cf) halfinteger. 

(b) In addition, there can be negatíve eigenvalues. To find them we suppose that 

the loop wavefunction is given by the first part of ( [3.1|) with k = in and the 
halfline part is pe~ KX , k > . The boundary conditions then yield a system 
of equations for f3, p which can be solved provided 

(cos <3> — cos kL) = a + . (3.8) 



sinh kL 1 + /i/t 

It is easy to see that under the condition a > \ (cos$ — 1) has no solution 
if fi > and a single root otherwise; in the case a < % (cos $ — 1) one more 
eigenvalue is added. 



4 A digression: a duality for graph scattering 

At this point we want to make a small detour to describe a general method to 
treat scattering problem on graphs. Recali that there is an equivalence between the 
spectral problem for one-dimensional Schròdinger operators with point interactions 
and certain Jacobi matrices which is known in the literature as a "French connection" 
0, 0, [| ^T], |3T|. We have been able to extend this duality recently to a wide 
class of Schròdinger operators on graphs JÏ4|]; here we want to illustrate that the 



same method is applicable to scattering problems. 



4.1 Schròdinger operators on a general graph 

Let us first collect so me notion we shall need to formulate the result; for more 
details we refer to |TJ]]. A graph F consists of a finite or countably infinite number 
of vertices V = {Xj : j G /} and links (edges) £ = {Cj n : (j,n) el/; Cix /}. 
Without loss of generality we may suppose that each pair of vertices is connected 



6 



by not more than one link; otherwise we just add some number of extra vertices. 
We assume that T is connected, so the set N(Xj) = {X n : n G v{j) C I \ {j}} of 
neíghbors of Xj , i.e., the vertices connected with Xj by a single link, is nonempty. 
Throughout we shall assume that J\f(Xj) is finite for any j G I . 

The graph boundary B is the subset of vertices having a single neighbor; it may 
be empty. We use the symbols Jg and Jj for the index subsets in I corresponding 
to B and the graph interior X := V\B , respectively. T has a /oca/ mètric structure 
coming from the fact that each link Cj n can be mapped to a line segment [0,^ n ] . 
It is also possible to equip the graph naturally with a global mètric by identifying it 
with a subset of a plane or a higher dimensional Euclidean space. The two mètrics 
may differ at a single link; the local one which is important for us is usually given 
by the arc length of the curve segment representing Lj n . 

Using the local mètric, we are able to introduce the state Hilbert space in 
the way we did it for the lasso graph and similar problems, namely as L 2 {T) := 
®(J,n)eij: -^ 2 (0> £jn) ■ Its elements, i.e., the wave functions, will be written as = 
{ipjn '■ (j,n) G Ic} or simply as {if)jn} ■ We shall suppose that the particle liv- 
ing on T is exposed to a potential; it is only important to know its vàlues on the 
graph links, i.e., a family of functions V := {Vj n } ; since we do not want deal with 
mathematical subtleties here, we suppose that all of them are essentially bounded, 
Vj n G L°°(0, £j n ) . Then we are able to define the operator H a = H(T, a, V) by 

H a {i> jn } := { + V jn ^ jn : (j, n) G I c } (4.1) 

with the domain consisting of all ip with ipj n G W 2 ' 2 (0, £j n ) subject to a set a of 
boundary conditions at the vertices which couple the boundary vàlues 

we have identified here x = with the vèrtex Xj . In general, there is vast family of 
boundary conditions which make the operator ( |4.1| ) self-adjoint. It can be character- 
ized by 4M 2 real parameters, where M is the number of graph links 0, [T^, ^] , and 
even if we restrict to local boundary conditions which do not couple the boundary 
vàlues belonging to different vertices, the number is still too large. 

As above we restrict ourselves to the simplest situation when the links connected 
in a vèrtex Xj satisfy the ó-coupling condition, i.e., ipj n (j) = 4>jm(j) ='■ i>j for all 
n, m G u(j) , and 

£ ^n(j) = '•,<••., (4-3) 

n£u(j) 

with a real-valued parameter atj G M (coupling constant). However, the results 
derived below can be reformulated easily for the case when (|4.3|) is replaced by a 



o"'-coupling or another type of local boundary conditions |Ï2|, [T] 

As in the particular case discussed in the previous sections the relation Q4.3| ) and 
other local couplings have an illustrative meaning of probability current conservation 
at the vèrtex; in a sense they represent an analogy of Kirchhoff 's law. This means, in 
particular, that they are independent of the lengths of the involved links. Moreover, 



7 



since the probability current is connected with the kinetic part of the Schrodinger 
equation, the coupling is also independent of the potential V as long as the latter 
is regular, which is the assumption we have adopted. At the graph boundary we 
employ the usual conditions 

if)j cos ujj + i/jj sin ujj = (4.4) 

with a parameter ujj ; integer and halfmteger múltiples of tt correspond to the 
Dirichlet and Neumann condition, respectively. 

4.2 Coupling two link bundles 

In the next step we attach a certain number of semiinfmite links to F which will 
support asymptotic solutions; in the Standard stationary picture we shall consider 
a combination of a falling and transmitted/reflected plane wave on each of them. 
We might regard these "external" links as a part of the graph boundary; however, 
it is convenient to treat them separately. A reason for that is the following: while 
we declared the intention to formulate the result for graphs with 5-couplings, it is 
desirable to have a coupling between the infernal and external links which is slightly 
more general than ((4.3|). This could be useful, e.g., if we want to study perturba- 
tively resonances which arise when eigenvalues of the original graph operator become 
embedded into the continuous spectrum of the leads. 

As another preliminary, therefore, consider two bundles of leads which support 
wavefunctions {f n }n=i an d {g m }m=i i the endpoints are placed to the point x — . 
Suppose first that we have separate 5-couplings for each bundle, 

/i(0) = ■ • ■ = /at(O) =: /(0) , 9l (0) = ■■■ = g M (0) =: g(0) , (4.5) 

together with J2n=i fn(®) = a f(®) an d Em=i 9n(®) = ^(O) ■ To couple the two 
bundles, we preserve the separate continuity ( |4.5|) and replace the derivative condi- 
tions by 

N M N M 

/(0) = «- 1 E fn(0) + 7 E ^(0) , 9(0) = 7 E /n(0) + à- 1 E ^(0) (4.6) 

n=l m=l n=l m=í 

with a complex parameter 7 ; an elementary integration by parts then shows that 
the corresponding boundary form vanishes under these conditions. The parameter 
modulus is the coupling strength; if the coupling is required to be time-reversal 
invariant, 7 has to be real. An overall í-coupling is achieved, of course, if at = 
à = 7 _1 . 

4.3 The S— matrix equation 

Suppose now that a bundle of rrij halflines, < rrij < 00 , is attached to the point 
Xj of T; the coupling being given by (|4.5|) , (|4.6|) with the parameters atj for the 
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graph links joined at Xj , àj for the external links at Xj , and 7j . We call the 
j'-th bundle £j and £j m will be the m-th halfline in it, so the full state Hilbert 
space will be now L 2 (T) © (®jei 0m=i L 2 (£jrn) • For the sake of brevity, the graph 
extended by the external links will be denoted as r e = T U £ , for the state Hilbert 
space we will use the shorthand L 2 (T e ) . The symbol H a = H(T e , {ctj, àj, jj}, V) 
means a Schròdinger operator on T e with the described coupling; for simplicity we 
assume that the potentials on the external links are zero. 

As usual the stationary scattering problem means finding a generalized eigen- 
vector of H a with prescribed behavior in the asymptotic region, i.e., a solution to 
the equation 

H a ^ = k 2 ^, (4.7) 

which belongs locally to D(H a ) satisfying all the domain requirement (in particular, 
the boundary conditions at each vèrtex) apart of the global square integrability, and 
such that 

1> jm {x) = a jm e- lkx + b jm e ikx (4.8) 

holds for x G £j m ■ The vectors a = {dj m } and b = {bj m } of dimension card£ = 
J2jei m j represent the incoming and outgoing amplitudes, respectively; we are in- 
terested in the operator that maps the former into the latter, b = Sa . 

To proceed further, we need some more notation. The symbol will denote 
the decoupled operator obtained from H a by changing the conditions (|4.3| ) at the 
points of graph interior X to Dirichlet, while at the boundary they are kept fixed; 
we also define K, a := {k : k 2 G a(H^), Imk > 0}. Next we take an arbitrary 
link C n j = [0, £jn) of T , the right endpoint being identified with Xj , and denote by 
"Ujni Vjn the normalized Dirichlet solutions to the corresponding component — /" + 
Vjnf — k 2 f of the Schròdinger equation ( [4.7[ ). In other words, we demand that the 
following boundary conditions are satisfied, 

uMn) = l-u' jn (£ jn ) = , v jn (0) = l-v' jn (0) = , (4.9) 

provided n G Ix ', at the graph boundary we replace the last requirement by Vj n (0) = 
sinu n and v'j n (0) = —cosu n . The Wronskian of these solutions equals 

W jn = -v jn (£ jn ) = u jn (0) (4.10) 

for n G Ix and Wj n = —Uj n (0) cos u n — u'j n (0) smu n otherwise. All these quantities 
depend in general on the spectral parameter k but we shall not indicate this fact 
explicitly. Now we can formulate the mentioned result: 

Proposition. Let k G" JC a with k 2 G M, Im k > . Under the assumptions given 
above, the corresponding on-shell scattering matrix for the graph T e is given by the 
following system of N := card / + card£ equations 

EV'n / V jn(^jn) \ . , 

TTr L· '\kt a ó Ü ~ ika jl jmjbj 1 

n£v(j)nl x VV J n \n&j(J) VV Í n ) 
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(4.11) 



—ikotj^j mjCXji — 2 



m=2 / 



íïí( y — - y 

\neu{j)ni x VV 3 n neu{j) vv i n 



à m Z_/ Ti7 22 tí7 m+ bjxiàj-ikrrij 



(4.12) 



,i (o? + ikrrij) — 2ik 



m=2 



and 



&jm = &ji + aji - a jm , m = 2,...,mj. (4-13) 



Remarks. (a) If N < oo the above relations represent a system of linear equations. 
In the opposite case they have to be interpreted as the appropriate operator equation 
on £ 2 . This can be done under some additional assumptions on T , e.g., if there are 



positive numbers ci, C2 such that c± < £j n < C2 holds for all (j, n) E Ic — see |]T4 
for more details. 

(b) The results generalizes easily to the situation when T as a subset of IR U 
is placed into a magnètic field, not necessarily homogeneous, described by a vec- 
tor potential A. The boundary conditions (|4.3|) are modified replacing if)j n (j) by 
ipj n (j) + iAj n (j) , where Aj n (J) is the tangent component of A to Cj n at Xj 
The particle abiding on T is supposed here to be an electron; otherwise A has to be 
replaced by —qA where q is the particle charge. The magnètic case can be handled 
by means of the unitary operator U : L 2 (T) — > £ 2 (r) which acts as 




(U^) ín (x) := exp i / A jn (y) dy fi 



X 



the vàlues fixed reference points. Then the functions (Uip)j n satisfy ( |4.3| ) 



and it is sufficient to replace the function vàlues ip n in (|4.11| ), (|4.12|) by e lAn ip n 
provided the magnètic phase factors Aj are chosen to obey the natural consistency 
condition 

Aj - A n = / A jn (y)dy. 

required by the wave function continuity. 

(c) Consider a simple situation when a single halfline is attached to every point 
of T and denote the "graph part" of the above system, i.e., the operator represented 
by the two sums at the Ihs of ( |4.11 ) as h . If the coupling is ideal, aj = for all 



j E I , the S-matrix is given by 

S h + ik 
h — ik 
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It is illustrative to compare this to the the formula used recently by Sadun and 
Avron |33| in a study of scattering on discrete graphs; the only difference is the 



replacement of — ik by e lh , the energy being 2cos/í; in this case. 

To prové the proposition, it is sufhcient to use the transfer matrices which re- 



late the Schròdinger equation solutions at both ends of each link ||14j| . Since the 
Wronskian is nonzero for k ^ K a , we get 

Í>i '■= ^jn(j) = u' jn (0)ïp n + v jn (£ jn )ífj' jn (n) , 

VVjn 

the sign change at the Ihs of the last condition reflects the fact that (|4.2|) defines the 
outward derivative at Xj . We express ^j n (n) from the first relation and substitute 
to the second one. This yields 



%n{3) w jn + w jn % 



for n G 1% , while at the graph boundary we get with the help of ( |4.4j ) instead 

vvjn 

Now one has just to substitute these vàlues into the boundary conditions at each 
vèrtex to arrive at the relations ( |4.11| )-( f4.13| ). | 



It is not difficult to check that the lasso graph with the 5-coupling can be treated 
within this general scheme. We use the normalized Dirichlet solutions at both loop 
"ends" , k~ l sin kx and — k^ 1 sin k(x—L) , and add a vèrtex into an "interior point" . 
Using ( |4.11| ) and ( |4.12| ) and excluding the function vàlues at the added point, we 
arrive after a straightforward calculation at the equations 

/ , e^ + e"** n , cos kL \ 

—k — — + Ik — + a \ip — iko = —%ka , 

\ sm kL sin kL J 

ip = b + a , 

from which we recover the reflection amplitude ( |3.4j) . 

5 Resonances 
5.1 The resolvent 

Let us return now to our model. The most natural way to study spectral properties 
of an operator is through its resolvent, and therefore we want to find it for H a (B) . 
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The "decoupled" resolvent is found easily: it is a matrix integral operator with the 
kernel 

/ e -iA(x-y) sin kx < sinfcQ>-L) g \ 

k sin kL 



Goo(a:,y; k) 



\ 







sinfca;< exp(ikx > ) 
k 



(5.1) 



where x< and x > mean conventionally the smaller and larger of the variables x, y , 
respectively. We abuse here again the notation and employ the same symbol for the 
arc-length variable on the loop and the lead as well as for the pair of them. 

Since Ha^u and iïoo are both self-adjoint extensions of the same symmetric 
operator with the deficiency indices (2, 2) , the resolvent of the former is by Krein's 
formula [Q, App.A] given by 

2 



G 



y;k) = G OQ (x,y;k) + £ X je F j (x)F^(y) , 



(5.2) 



where the symbol "t" means transposition, Fj are vectors of the corresponding 
deficiency subspaces which we shall choose in the form 



with 



F 1 (x) :-- 



w\x) 



w{x) 
. 



Fo(x) :-- 







(5.3) 



Ax e sin kx — sin k(x — L) 
sin kL 



and Xji are coefficients to be found. Introducing 

pL roo 

hi := / w(y)v(y)dy } h 2 := / e iky g(y)dy 
Jo Jo 

for a given (^j £ TC , we find easily that the boundary vàlues of the function (^j :- 
(H a) ^ u - k 2 )" 1 Q are in view of (|5^ ) given by 

m(0) = w(L) = An/ii+ Ai 2 /í2, /(0) = À 2 i/ii+ A22/Í2, 



«'(O)-u'(L) = h x + 



2fc 



cos$ — cos kL) (Àn/ii+ \12h2) 



sin fc-L 

f(0) = / l2 + z fc(A 21 / il +A 22 / Í2 ) . 

However, Cfj belongs to D(H a fJ-LJ ) for any (^j , so substituting these boundary 
vàlues into [\2.7\) we get a system of four linear equations which yields the sought 
coefficients: 

1 — ifik uj 
— p— . Al2 = "p' 

(5.4) 



A 



11 



A 



21 



V 



X 



l·l 



2k 



22 



cosQ— cos kL 
sin kL 



a 



Lü 
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with 

[ cos<£> — cos kL 

V = Via, uw; k) := (1 — iiik) 2k — —a 

sm kL 



icü 2 k. (5.5) 



In the case of 5-coupling, // = 0, uj — 1, the coefficients acquire a particular ly 
simple form, Xj£ = — Z> -1 , j, í = 1, 2 . 

5.2 Pole trajectòries 



As usual in such situations [|, [Ï8[ the singularities of Goo(x, y; k) cancel with those of 
the added term in fl5.2p and the resolvent poles are given by zeros of the denominator 
(|5.5| ); the exception is represented by the case of an integer or halfmteger <fi . 

For the sake of simplicity, we shall speak mostly about the 5-coupling situation. 
If the coupling is ideal, a = , the pole condition becomes 

2(cos $ — cos kL) = —i sin kL (5-6) 

and one is able to solve it explicitly. No singularities exist in the upper halfplane, 
hence we write 

k = k — in . (5.7) 

Substituting into the above condition, we find that for |$| < | (mod ir) there is a 
pair of poles with k = irn/L and 
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- In (2(-l) n cos $ ± V4 cos 2 $ - 3) , (5.É 



where (— l) n cos $ > . On the other hand, for the remaining vàlues of $ the poles 

ln 3 
2L 



are found at the line parallel to the real axis with V = ~^tr an d 




k — ± — ■ arccos —= cos$ . (5.9) 



We see that both poles are in the open lower halfplane with the exception of $ = nu , 
í.e., (f) integer or halfmteger, when one of them turns into an embedded-eigenvalue 
pole at the real axis. The pole trajectòries with respect to $ are not smooth 
despite the analytic form of the condition ( |5.6| ); this is due to the fact that V = 
at the crossing points ^ (rm — | ln 3J , so the implicit-function theorem does not 
apply there. A similar picture is obtained for the boundary conditions ( |2.7| ) with 
a — fj, — and \to\ < y/2, in which case the "horizontal" line has r\ — | ln . On 

the other hand, in the case |a>| > y/2 the pole trajectòries are "vertical" segments 
with k = im/L only. 

If a 7^ the 5-coupling pole condition (|5.6|) is replaced by 



2/c(cos $ — cos kL) — (a — ik) sin kL . 
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Figure 2: Pole trajectòries from the condition (5.10) for different vàlues of the 
coupling constant (dashed: a = 0.5 , full: a = 0.1 , dotted: a = 0.05 ) 



Writing separately the real and imaginary parts with the help of the parametrization 
(|5.7|) , we find that for k — irn/L a zero can exist only at the real axis if $ = mn . 
For other vàlues of k the pole condition can be cast into the form 

277 — k cot kL . ^ . 

cothriL = 2+ a— ' r -, (5.10) 

T][r] — a) + k 

which has to be solved numerically. The resulting pole trajectòries are shown on 
Figure 2. 



6 Decay of loop states 

Up to now we have considered the lasso graph as a scattering system. Now we shall 
suppose that the system is prepared at an initial instant in a state the wavefunction 
of which is localized at the loop. It is not so important how such a situation is 
realized. For instance, one can place an electron at an isolated ring and "switch in" 
the junction at t = . The state is generally unstable under the evolution governed 
by and we are interested in the way in which it decays. 

Since we have an explicit expression for the resolvent, we are able in principle to 
write the non-decay amplitude explicitly [JÏT], Sec. 3.1]. However, instead of trying 
to evaluate this function we limit ourselves to elucidation of its bàsic properties. 



6.1 Spectral decomposition 

The relations ( |5.2| ) and ( |5.4j ) imply, in particular, that the resolvent form z \— > 
(^,(H atliiU) — z)" 1 is a meromorphic function including its continuation to the 
second sheet. Its possible poles are associated with the discrete spectrum of H a41UJ 
which we also know explicitly. Since these are the only singularities, the function 
(ifj, (H a l·i w — -) _1 ipj is analytic for all ip belonging to the complement / Hp{H a41LÚ ) L . 
In particular, it is uniformly bounded in any finite part of the strip |Imz| < 1 , and 
thus by the bàsic criterion |!2], Thm.XIII.19] such a vector belongs to Ti. ac (H aíl·liüJ ) . 

Consequently, our Hamiltonian has no singularly continuous spectrum. The 
initial state can be therefore decomposed into ip = ip p + ip ac and the corresponding 
non-decay amplitude equals 

(i>, U t tfj) = (V> P , U t ^j p ) + (V»oc, U t ^ ac ) , (6.1) 

where U t := exp{— iH a41Lü t} . The second term on the rhs goes to zero as t — > oc 
in view of the Riemann-Lebesgue lemma; the first one is a linear combination of 
exponentials with coeficients coming from the Fourier decomposition of ip p . If just 
one of them is nonzero, then the decay law of the corresponding loop state given 
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by squared modulus of ( |6.1[) has a finite nonzero limit. Such a behavior is typical 
for unstable systems having bound states with a nonzero Fourier component in a 
decaying state; it has been observed recently in another context — see [fL5|| , and also 



|Ï9| where, however, the effect may be also related to the threshold violation of the 



Fermi golden rule discovered by Howland [29]. 

If the loop state contains a superposition of a larger number of eigenvectors, the 
nondecay probability does not go to zero as t — > oo but a limit does not exist. In 
view of the above discussion, such a situation can occur in the present model only if 

(a) there are two negative eigenvalues (see Remark (b) at the end of Section 3), or 

(b) if $ = nir with n e Z . The asymptotic behavior of the decay law depends then 
on the coupling parameters. If all the involved eigenvalues are commensurate, the 
asymptotics is periòdic; this happens always if there is no negative-energy bound 
state. In the general case the decay law asymptotics is quasiperiodic. 



6.2 What has all this in common with neutral kaons? 

Concluding this study, let me mention one more tòpic to which Larry Horwitz made 
a contribution, namely the decay theory of neutral kaons. This subject attracted 
attention at the end of the sixties as an example of a system with a substantially 
nonexponential decay law exhibiting different time scales, as well as the possibility 
to "recreate" decayed partides by performing a set of noncompatible measurements. 

Mesoscopic physics makes it possible to tailor systems in which similar effect can 
be observed. Consider our lasso graph with the initial wavefunction u on the loop 
such that x >— > e lAx u(x) has no defmite symmetry with respect to the connection 
point x = (say, u(x) = e - tx ( A - 27rn / L ) y if the flux value is integer, the A— 
even component represents a superposition of embedded-eigenvalue bound states 
and thus it survives, while the A-odd one dies out. In a real life experiment, of 
course, we cannot ensure that is exactly an integer, hence we shall have rather 
a fast and a slowly decaying part of the wavefunction; recali the pole trajectòries 
discussed in Section 5.2. 

Moreover, consider a loop to which two halfline leads are attached at different 
points and assume that we are able to switch the coupling in and out independent ly. 
We wait until the A-oàà part in the above described experiment essentially decays 
while the longliving component is still preserved, and switch from the first lead to 
the second one. Now the symmetry with respect to the other junction is important. 
If the surviving part of the wavefunction is a superposition of an A-even and an 
A-odd part with respect to the latter, the scenario repeats. Of course, the "second 
decay" may produce a smaller component A-odd with respect to the first junction, 
so the analogy is complete. 
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Figure captions 

Figure 1. A lasso graph 

Figure 2. Pole trajectòries from the condition (5.10) for different vàlues of the 
coupling constant (dashed: a — 0.5 , full: a — 0.1 , dotted: a = 0.05 ) 
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